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It is pointed out that there exists an interesting strong and weak duality in the Landau-Zener-Stueckelberg 
potential curve crossing. A reliable perturbation theory can thus be formulated in the both limits of weak and 
strong interactions. It is shown that main characteristics of the potential crossing phenomena such as the Landau- 
Zener formula including its numerical coefficient are well-described by simple (time-independent) perturbation 
theory without referring to Stokes phenomena. A kink-like topological object appears in the "magnetic" picture, 
which is responsible for the absence of the coupling constant in the prefactor of the Landau-Zener formula. It is 
also shown that quantum coherence in a double well potential is generally suppressed by the effect of potential 
curve crossing, which is analogous to the effect of Ohmic dissipation on quantum coherence. 



1. Introduction 

The potential curve crossing is related to a wide 
range of physical and chemical processes, and the 
celebrated Landau-Zener formula correctly 
describes the qualitative features of those pro- 
cesses P|-|lq| . It has been recently shown [O] that 
the potential crossing problem contains interest- 
ing modern field theoretical ideas, namely, the 
duality and gauge transformation. 

The adiabatic and diabatic pictures in poten- 
tial curve crossing problem are related to each 
other by a field dependent su(2) gauge transfor- 
mation and we point out that this trans- 
formation leads to an interchange of strong and 
weak potential curve crossing interactions, which 
is analogous to the electric and magnetic duality 
in conventional gauge theory p2| . This strong 
and weak duality allows a reliable perturbative 
treatment of potential curve crossing phenomena 
at the both limits of very weak (adiabatic pic- 
ture) and very strong (diabatic picture) potential 
crossing interactions. 



2. A model Hamiltonian of potential curve 
crossing and duality 

To analyze the potential curve crossing, we 
start with a model Hamiltonian defined in the 
so-called diabatic picture p,p| 
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P 



V 1 (x)+V 2 (x) , V 1 {x)-V 2 (x) 



-cr 3 



(1) 



where 03 and o\ stand for the Pauli matrices. We 
assume throughout this article that the potential 
crossing occurs at the origin, Vi(0) = V^O) = 
(see Fig. 1). 

If one neglects the last term in the above Hamil- 
tonian, one obtains the unperturbed Hamiltonian 
in the diabatic picture 



1 

2m 



V 1 (x) + V 2 (x) , V x {x)-V 2 {x) 



03.(2) 



* Invited talk presented at 7th Asia Pacific Physics Confer- 
ence, August 19-23, 1997, Beijing, China (to be published 
in the Proceedings) 



This Hamiltonian Hq describes two potentials, 
which are decoupled from each other. The last 
term in ([!]), Hj = o~\/g with a constant g, causes 
the transition between these two otherwise in- 
dependent potential curves. In other words, if 
one takes g — > large, this case physically cor- 
responds to a complete potential crossing from 
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Figure 1. Landau-Zener process in the diabatic 
(electric) picture. 



Figure 2. Landau-Zener process in the adiabatic 
(magnetic) picture. 



a view point of adiabatic two-potential crossing 
in Fig. 2. Namely, g stands for the strength of 
potential crossing interaction, and g — > large cor- 
responds to a very strong potential crossing inter- 
action. On the other hand, if one lets g smaller, 
the effects of the last term in (1) become substan- 
tial and the Hamiltonian Hq ( 2 ) does not present 
a sensible zeroth order Hamiltonian. 

To deal with the case of a small <?, we perform 
the non-Abelian "gauge transformation," 



H l _ e i6{x)a 2 /2 He -ie(x)a 2 /2 



(3) 



where a 2 is a Pauli matrix. The Hamiltonian in 
the new picture is given by 



2m 



P - -^d x 8{x)a 2 



V 1 (x) + V 2 (x) 



Vi(x)-V 2 (x) 1 . . 

cos 0(x) H — sin v(x) 

2 9 



^4) 



V 1 {x)-V 2 (x) 



smf (x) 



■ cosv(x) 



f) 



a i. 



To eliminate the potential curve mixing, the last 
term of (|5|), we choose the gauge parameter 8{x) 

as a 



C ote( x ) = 9 v ^- v ^ =m. 



(5) 



We then obtain the Hamiltonian in the adiabatic 
picture 

H' = H' + Hi, (6) 
where 



H' 



and 



-P 



U 1 {x) + U 2 {x) 



2m 2 
U 1 (x)-U 2 (x) 

o CT 3, 



(7) 



4m 



\pd x 6(x) + d x 9(x)p) a 2 



(8) 
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The potential energies in the adiabatic picture are 
related to those in the diabatic picture as (Fig. 2) 



V 1 (x)+V 2 (x) 



± 



V 1 (x)-V 2 (x) 



(9) 



From the definition of the gauge parameter in (||) , 
the "gauge field" d x 9(x) is expressed as 



d x e{x) 



(10) 



The transition from the diabatic picture to the 
adiabatic picture is a field dependent transforma- 
tion. 

In the adiabatic picture, the o~ 2 dependent term 
in the interaction H'j (^) causes the potential 
crossing. If one neglects H'j, the two potentials 
characterized by U\{x) and U 2 {x) do not mix 
with each other: Physically, this means no po- 
tential crossing. This suggests that H'j is propor- 
tional to the coupling constant g, since a small g 
corresponds to weak potential crossing by defini- 
tion. This is in fact the case as is clear from (10) 
and (|). 

We thus conclude that the two extreme limits 
of potential crossing interaction should be reli- 
ably handled in perturbation theory; namely, the 
strong potential crossing interaction in the dia- 
batic picture, and the weak potential crossing in- 
teraction in the gauge transformed adiabatic pic- 
ture. This is analogous to the electric-magnetic 
duality in conventional gauge theory Q : The di- 
abatic picture may correspond to the electric pic- 
ture with a coupling constant e = l/g, and the 
adiabatic picture to the magnetic picture with a 
coupling constant g. 

A general criterion for the validity of perturba- 
tion theory in the adiabatic picture (o) is 



-\d x 9(x)\ « \p(x)\, 



(11) 



which is expected to be satisfied when the cou- 
pling constant g is small and the incident particle 
is sufficiently energetic. 



3. Landau-Zener formula 

As an illustration of the duality discussed in 
Section |^, we re-examine a perturbative deriva- 
tion of the Landau-Zener formula in both of the 
adiabatic and diabatic pictures For def- 

initeness, we shall assume V((0) > V 2 (0) as 
in Fig. 1. 

Let us start with the adiabatic picture with 
weak potential crossing interaction. Since the 
gauge field generally vanishes, d x 9(x) — > for 
| a; | — * oo, we can define the asymptotic states in 
terms of the eigenstates of H' (^) . We define the 
initial and final states $i and <&/ by 



*i(a:) 



<px{x) 






<p 2 (x) 



(12) 



which satisfy 



±f + U l{ x) 
±f + U 2 (x) 



<p 2 (x) 



Etpi(x), 
Eip 2 (x). 



(13) 



We then obtain the potential curve crossing prob- 
ability due to the perturbation H'j (0) 



?7T 

w(i^f) = -\($ f \H' I \$ i )\ 2 . 

The transition matrix element is given by 

h 2 



(14) 



(*f\H'i\*i 



4m 



dx d x 9(x) 



(15) 



x [-^' 2 {x)^Pi{x) + ^(x^ix)} 



To evaluate the matrix element, we use the 
WKB wave functions B: 



cos 



exp 



1 f x 

-- J dx\p-L(x)\ 

dxpi(x) - - 



\/pi(x) 

for x > ai and x < ax, respectively, and 
C 2 



(16) 



, 2 cc 



exp 
1 

Ti 



l r 



dxp 2 (x) 



(17) 
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for x > a 2 and x < a 2 , respectively. The semi- 
classical momenta in the adiabatic picture are de- 
fined by 



pi, 2 (x) = ^2m[E-U h2 (x)}, 



(18) 



and a\ and a 2 denote the classical turning points 
(see Fig. 2). The normalization of <pi(x) is chosen 
as C\ = 2^/m to make the probability flux of the 
incident wave unity. On the other hand, the final 
state wave function in ( |l4| ) has to be normalized 
by the delta function with respect to the energy 
($2 1 $2) = 5{E' 2 - E 2 ) and this specifies C 2 = 
2^/V2nh. 

We estimate the matrix element (15) by using 
the oscillating parts of the wave functions (|l6| ) 
and (|l7|). This treatment is justified if the follow- 
ing conditions are satisfied: 

(i) |p(0)| — > large and m — > large with u = 
\p(0)\/m kept fixed such that non-relativistic 
treatment is valid in the physically relevant re- 
gion. 

(ii) g — > small, but with 



-«^b(0)l 2 

g Am 



to ensure (11) and the condition 



(3 < a, 



(19) 



(20) 



where a is an average turning point and is a 
typical geometrical extension of d x 9(x). If (|2(]) is 
satisfied, we can estimate the matrix element by 
using the oscillating parts of wave functions only 
since d x 8(x) rapidly goes to zero for |x| ^> (3 on 
the real axis. 

The integral (15) is then written as 



ihC\C 2 
8m 



dxd x 6(x) (21) 



cxp 



exp 



f x i f x 

I dxpi(x) — — I dxp 2 (x) 

J ax "> Ja 2 

dxp 2 (x) 



dxpi(x) + — 

n 



where we have set pi(x)/p 2 (x) = 1 in the prefac- 
tors. This is justified if U/\p(0)\ <C /3, the char- 
acteristic length scale of the present problem, by 



letting p(0) large as is specified in (i) . Therefore 
we need to evaluate an integral of the form 



I = 



dx d x 9(x) 



(22) 



x exp 



- / dxpi(x) - - dxp 2 

J at n Ja 2 



We here present an explicit evaluation of (22) 
for the linear potential crossing problem, Vi(x) = 
V[(0)x and V 2 (x) — ^'(O)^, on the basis of lo- 
cal data without referring to Stokes phenomena. 
For sufficiently large energy, E — {U\ + U 2 )/2^> 
(Ui — U 2 )/2, the difference of momenta can be 
approximated as [see @], 



dx[pi(x) -p 2 {x)] 



(23) 



dx ■ 



lo v(x)g(3 
where we used 



V^ 2 + P 2 



v(x) 



9 — 2 — x = p> 



(24) 



1 



'2m 



E 



U^ + U^x) 



and v(x) is approximated to be a constant v 
v(0) in the following. We also have from ( [To|) 



d x 9(x) = - 



x 2 + (3 2 



(25) 



and thus 

I ~ — exp 

dx 



if if 

-J dxp 1 (x)--J dxp 2 (x) 







exp 



x 2 + (3 2 




exp 



2i 

1ivgl3 J 

r0 



dx yjx 2 + /3 2 

^ dxpi{x) - ^ / dx p 2 (x) 
h Jat h J a2 



/oo 
dx exp[— iaF(x) 
-00 



where 



Fix) = / dx \Jx 2 + 1 + — 
Jo l0L 



ln(x 2 + l), 



'0 

2/3 4 
Jw~g ~ hvg 2 [V{(0) - Vfipj] 



>0. (27) 



(26) 
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F'{x) = Vx 2 + 1 + — 



We evaluate the integral (|26|) by a saddle point 
approximation with respect to a. We thus seek 
the saddle point 

-J^T = 0, (28) 

which is located between the real axis and the 
pole positions x — ±i of d x 9(j3x) so that we can 
smoothly deform the integration contour; these 
poles also coincide with the complex potential 
crossing points. If one sets x = iy in ( p8|) for 
— 1 < y < 1, one has 

1 2y 



which has a unique solution 



i 2 
i + - [- 
2\a 



2/3 



(29) 



(30) 



for large a. (The complex conjugate of x s is lo- 
cated in the second Riemann sheet.) For this 
value of the saddle point 

2/3 



F(x s ) = 



F"(x s ) 



dx \fx? 



1 



7ri 2 i 1,(2 
1 1 In - 

4 3 a la. \ct 



2/3 



-3i 



(!) 



1/3 



(31) 



We thus have a Gaussian integral which decreases 
in the direction parallel to the real axis 



x exp 



~i f if 

- / dxp^x)-- I dxp 2 (x) 

-3(f) 4/ V* s ) 2 

/ 7ra\ 
V 4/ 

- / dxpi(x)-^- / dxp 2 (x) 

ft Ja, ft Ja 2 



dx exp 



x exp 
From (E2l) we obtain 



(®fW'i\®i 



,2/3 



x exp 



M 

dxpi(x)— / dxp 2 (x) 

til •» a 2 

7T 



^ 2 [y/(o)-y 2 '(o)] 



(33) 



It is interesting that the numerical value of the 
coefficient of the above expression, y^e 2 / 3 /\/3 = 
1.99317, is very close to the canonical value 2 
and we replace it by 2 in the following. As for 
the past analysis of the prefactor in the time- 
dependent perturbation theory, see papers in J13| . 
We thus have the transition probability from (111) 



w(i -> /) 



4 sin 



x exp 



dxpi(x) 



2tt 



dxp 2 (x) 



2 exp 



hvg 2 [vm-vm] 

2vr 

~hv 9 2 [vm-v>m 



(34) 



where we replaced the square of sine function by 
its average 1 /2 in the final expression. We empha- 
size that the numerical coefficient of w(i — > /) 
is fixed by time-independent perturbation the- 
ory and the local data without referring to global 
Stokes phenomena; this is satisfactory since linear 
potential crossing is a locally valid idealization. 

We interpret that w(i — > /) in (34) expresses 
twice the non-adiabatic transition probability. 
Notice that our initial state wave function con- 
tains the reflection wave as well as the incident 
wave. Therefore the transition probability per 
one crossing is given by the half of (p4|) , 



exp 



2- 



hvg2[V{(0)-Vl(0)] 



(35) 



which is the celebrated Landau-Zener formula (J, 
D . Our perturbative derivation presented here is 
conceptually much simpler than the past works ||, 
mil, and it should be useful for a pedagogical 
purpose also. 

It is interesting to study the same problem in 
the diabatic picture in Fig. 1 with Hi = o\jg 
for large g. The evaluation of the matrix element 
is the standard one described in the textbook of 
Landau and Lifshitz B, for example. We have 



G 



for E > 0. 



>(i - /) 



^(q)5 2 [^/(o) - vm 

dxp 2 (x)-- / <ia;pi(a;)-- 



o 

4tt 



^ 2 [^i'(0)-F 2 '(0)]' 



4 
(36) 



[u(0) is the velocity at the crossing point, v(0) = 
^2E/m.] We again interpret (|3^) as twice the 
potential crossing probability because our initial 
state wave function contains the reflection wave 
as well as the incident wave. The transition prob- 
ability per one potential crossing is given by the 
half of ©. 

A simple interpolating formula, which repro- 
duces d34| ) in the weak coupling limit and ([36]) in 
the strong coupling limit, is given by 



w(i — > /) ~ 2 exp 



1 — exp 



2tt 



hvg 2 [V{(0)-Vl(0)] 
2tt 



hvg*\V{(0) - V a '(0)] 



(37) 



This expression is also consistent with the (semi- 
classical) conservation of probability 

Motivated by duality, we re-examined a pertur- 
bative derivation of the Landau-Zener formula, 
and we re-derived the formula (|34|) including its 
numerical coefficient on the basis of perturbation 
theory. However, our final result ( |34| ) in the adi- 
abatic picture does not contain the coupling con- 
stant as a prefactor. This is related to an interest- 
ing topological object in the present formulation. 



From the definition of (10), the "gauge field" sat- 
isfies the relation 



dxd x 0(x) = 0(oo)-0(-oo) 



(38) 



which is independent of the coupling constant g; 
we assume f{x) — > ±oo for x — » ±oo, respec- 
tively. Because of this kink-like topological be- 
havior of 0(x), the coupling constant does not 
appear as a prefactor of the matrix element in 
perturbation theory if the wave functions spread 
over the range which well covers the geometrical 



size of d x 6(x). The precise criterion of the valid- 
ity of perturbation theory is thus given by ([111): 
This condition is in fact satisfied if the conditions 
(|l9|)-(p0|) are satisfied. For small values of x, the 
small coupling g helps to satisfy (|ll|). Even for 
the values of x near the average turning point a, 
we have 

h dx e(a)\^U?-)*«-c\p(a)\, (39) 
2 2 \a ) a a 

where (3 stands for the typical geometrical size 
of d x 8(x). The estimate in the left hand side 
is based on linear potentials (p5|), but we expect 
that the condition is satisfied for more general po- 
tentials as well. We thus clarified the basic mech- 
anism why the prefactor of the Landau-Zcncr for- 
mula d35| ) should come out to be very close to 
unity in time-independent perturbation theory. 

4. Discussions 

Motivated by the presence of interesting weak 
and strong duality in the model Hamiltonian ([!]) 
of potential curve crossing, we re-examined a per- 
turbative approach to potential crossing phenom- 
ena. We have shown that straightforward time- 
independent perturbation theory combined with 
the zeroth order WKB wave functions provides 
a reliable description of general potential cross- 
ing phenomena. Our analysis is based on the 
local data as much as possible without referring 
to global Stokes phenomena. Formulated in this 
manner, perturbation theory becomes more flex- 
ible to cover a wide range of problems. 

The effects of dissipative interactions on macro- 
scopic quantum tunneling have been extensively 



analyzed in the path integral formalism |14| and 
also in the canonical (field theoretical) formal- 
ism [H| . It is generally accepted that the Ohmic 
dissipation suppresses the macroscopic quantum 
coherence; in fact, an attractive idea of a dissipa- 
tive phase transition has been suggested . 

It is plausible that the effects of potential 
curve crossing with nearby potentials influence 
the quantum coherence of the two degenerate 
ground states. One can in fact confirm that 
the potential curve crossing generally suppress 
the quantum coherence by using the perturbation 
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theory for both limits of strong and weak curve 
crossing interactions jll]]. From a view point of 
symmetry, the lowest order perturbation in the 
present problem and the dissipative interaction in 
the Caldeira-Legget model both correspond 
to a dipole approximation. However, a perturba- 
tive analysis of basically non-perturbative tunnel- 
ing phenomena requires a great care. In Ref. 
an explicit calculation of rather limited configura- 
tions has been performed, which in fact indicates 
the general suppression of quantum coherence by 
potential curve crossing. This suppression phe- 
nomenon of quantum coherence may become im- 
portant in the future when one takes the effects 
of the environment into account in the analysis of 
potential curve crossing processes. 

From a view point of general gauge theory, it is 
not unlikely that the electric-magnetic duality in 
conventional gauge theory is also related to some 
generalized form of potential crossing in the so- 
called moduli space jL2j. We hope that our work 
may turn out to be relevant from this view point 
also. 

REFERENCES 

1. L. D. Landau, Phys. Z. Sow. 1, 89 (1932); 2, 
46 (1932). 

2. C. Zener, Proc. Roy. Soc. A137, 696 (1932). 

3. E. C. G. Stueckelberg, Helv. Phys. Acta. 5, 
369 (1932). 

4. L. D. Landau and E. M. Lifshitz, Quantum 
Mechanics, translated by J. B. Sykes and J. 
S. Bell (Pergamon Press, 1977). 

5. M. S. Child, Molecular Collision Theory 
(Academic Press, 1974), and references 
therein. 

6. D. S. F. Crothers, Adv. Phys. 20, 406 (1971). 

7. J. B. Delos and W. R. Thorson, Phys. Rev. 
A6, 728 (1972); A9, 1026 (1974) (E). 

8. E. E. Nikitin and S. Ya. Umanskii, Theory 
of Slow Atomic Collisions ( Springer- Verlag, 
1984), and references therein. 

9. M. S. Child, Mol. Phys. 16 (1969) 313. 

F. C. Carter, Molecular Electronic Devices 
(Marcel Dekker, 1982). 

H. Nakamura, J. Chem. Phys. 97, 256 (1992). 

10. C. Zhu and H. Nakamura, J. Math. Phys. 33, 



2697 (1992); J. Chem. Phys. 97, 1892 (1992) 
102, 7448 (1995). 



11. K. Fujikawa and H. Suzuki, [hcp-th/9704062 
(to be published in Phys. Rev. A). 

12. P. A. M. Dirac, Proc. Roy. Soc. A133, 60 
(1931). 

C. Montonen and D. Olive, Phys. Lett. B72, 
117 (1977). 

N. Seiberg and E. Witten, Nucl. Phys. B426, 
19 (1994). 

13. A. M. Dykhne, Sov. Phys. JETP 14, 941 
(1962). 

M. Matsuzawa, J. Phys. Soc. Jpn 23, 1383 
(1967). 

J. P. Davis and P. Pechukas, J. Chem. Phys. 
64, 3129 (1976). 

14. A. O. Caldeira and A. J. Legget, Phys. Rev. 
Lett. 46, 211 (1981); Ann. Phys. (N.Y.), 149, 
374 (1983). 

A. J. Bray and M. A. Moore, Phys. Rev. Lett. 
49, 1545 (1982). 

S. Chakravarty, Phys. Rev. Lett. 49, 681 
(1982). 

G. Schoen and A. D. Zaikin, Phys. Rep. 198, 
237 (1990), and references therein. 

15. K. Fujikawa, S. Iso, M. Sasaki and H. Suzuki, 
Phys. Rev. Lett. 68, 1093 (1992); Phys. Rev. 
B46, 10295 (1992). 



